Magnetic orders and topological phases from f-d exchange in pyrochlore iridates 
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We study theoretically the effects of f-d magnetic exchange interaction in the R2lr20r pyrochlore iridates. 
The R 3+ /-electrons form localized magnetic doublets due to the crystal field environment, while the Ir 4+ d- 
electrons are more itinerant and feel a strong spin-orbit coupling. We construct and analyze a minimal model 
capturing this physics, treating the Ir subsystem using a Hubbard-type model. First neglecting the Hubbard 
interaction, we find Weyl semi-metal and Axion insulator phases induced by the f-d exchange. Next, we find 
that f-d exchange can cooperate with the Hubbard interaction to stabilize the Weyl semi-metal over a larger 
region of parameter space than when it is induced by d-electron correlations alone. Applications to experiments 
are discussed. 
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The discovery of time-reversal invariant topological band 
insulators has opened new terrain in the study of topological 
states of matter[ 1 —3 1 . As spin-orbit coupling (SOC) is essen- 
tial to realize some topological states, and given the propensity 
of electron correlations to induce a variety of novel phenom- 
ena, there is now intense interest in hd transition metal ox- 
ides with simultaneously strong (SOC) and intermediate elec- 
tron correlation. The pyrochlore iridates R 2 Ir 2 07, where R 
is a lanthanide (or Y), have been proposed to host various 
topological phases. Topological band insulators, topologi- 
cal Mott insulators [4], and Weyl semi-metal (WSM) statesj5], 
have been proposed||4]-|6]. In the WSM, which may occur as 
a magnetically ordered state in pyrochlore iridates(5] |6), the 
Fermi surface is composed of discrete, topologically stable 
Weyl points. In the related spinel compounds A2OS4O6, an 
Axion insulator (AI) ||5] 171- tTOl with bulk magnetic order and 
large magnetoelectric effect, has been proposed! 1 X | . 

The current theoretical analyses of pyrochlore iridates have 
focused either on the Ir subsystem|4-6 12 1, or (excluding 
non-magnetic R = Y, Eu) on the magnetic moments formed 
by localized R 3+ /-electrons fl3lfi7l . However, as we argue 
below, coupling between the local moments and the Ir 4+ d- 
electrons may be important in some compounds, but remains 
largely unexplored (see [18 19] for related prior work). More 
generally, most work on correlated, strong SOC materials is 
focused on correlations among e?-electrons. An alternate route 
to introduce correlation into a strong SOC system is to couple 
strong-SOC itinerant carriers to local magnetic moments, and 
such systems are certainly deserving of greater attention. 

In this Letter, we consider the effect of f-d exchange in 
pyrochlore iridates. The f-d exchange and the Ir electrons 
together generate a RKKY exchange between the localized 
moments at the R sites, which induces magnetic ordering on 
the R subsystem. The magnetic order on the R subsystem 
further modifies the electronic structure on the Ir subsystem, 
leading to WSM and Axion insulator phases. Notably, when 
f-d exchange is combined with correlation of the (^-electrons, 
we find the WSM is stabilized over a much wider region of 
parameter space than was found for d-electron correlations 
alone||6)- Based in part on this observation, we propose that 



Nd2lr207 is a candidate to realize the WSM. 

In many pyrochlore iridates, a metal-insulator transition 
and/or magnetic order occurs at a temperature scale of ~ 
100 K. Such a scale is presumably too large to be driven by 
RKKY interaction, and instead is probably set by d-electron 
magnetic exchange interaction, so the emphasis on Ir subsys- 
tem may be justified. However, in some compounds, inter- 
esting transport and magnetic properties only occur at lower 
temperatures. For instance, in Pr 2 Ir 2 07, a chiral spin liq- 
uid phase has been proposed at 0.3 K < T < 1.5 K to 
account for the anomalous hall effect observed in this tem- 
perature window, where no clear signature of magnetic order 
is observed[20]. In Nd 2 Ir 2 07, the metal-insulator transition, 
which seems to be associated with magnetic order, occurs at 
36K and can be suppressed by the application of pressure lEIT - 
l23l . It should be noted that other studies of the Nd compound 
show different behavior, apparently due to differences in sam- 
ple preparation [24 , 23 • 

We now describe our theoretical model for R 2 Ir 2 07. The 
R and Ir sites each form a pyrochlore lattice of corner-sharing 
tetrahedra, which can be viewed as a FCC lattice with four- 
site basis. The non-Kramers R 3+ ions (R = Pr, Tb, Ho) have 
an even number of /-electrons, while the Kramers ions (R 
= Nd, Sm, Gd, Dy, Yb) have an odd number of /-electrons. 
The large spin-orbit coupling of /-electrons then leads to a lo- 
cal magnetic moment with integer (non-Kramers) or half-odd- 
integer (Kramers) total angular momentum J. The (2 J + 1)- 
fold degeneracy is then split by the crystal field at the 
R-site. In the Kramers case, this leads to a magnetic doublet 
ground state. In the non-Kramers case, the crystal field splits 
the angular momentum multiplet into magnetic doublets and 
non-magnetic singlets. A magnetic doublet is the ground state 
for the R = Pr iridate lfl4l[T5l . as well as for isostructural insu- 
lating compounds such as Tb 2 Ti 2 07 and Ho 2 Ti 2 07||26] |27), 
so we assume a magnetic doublet ground state. We ignore 
effects of higher crystal field levels, which is valid for suffi- 
ciently large energy gap between the crystal field ground state 
and first excited state. This is reasonable for the R = Pr, Nd 
cases of greatest interest, where the reported gaps are 168 K 
and 300 K, respectively [ 28 29]. The R-site moment is thus 
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described in all cases by an effective spin- 1/2 pseudospin r. 

The d-electrons of the Ir subsystem are more itinerant. Ig- 
noring trigonal distortion of the oxygen octahedra surround- 
ing Ir sites, the 5d-electrons of Ir 4+ occupy the lower i 2ff man- 
ifold. The atomic SOC splits the 6-fold degenerate t 2g spin- 
orbital states into an upper doublet with j^f =1/2 and lower 
quadruplet with j e s = 3/2. In the tight-binding description 
of the electron tunnelings, a sufficiently strong SOC splits the 
Jeff = 3/2 bands from j e s = 1/2 so that j e s = 3/2 bands are 
completely filled and the j e g = 1/2 bands are half-filledH). 
Therefore, the Ir subsystem can be approximately viewed as a 
pyrochlore lattice system with one j & s = 1/2 doublet electron 
per Ir site. 

For the exchange coupling between R pseudospin r and 
Ir effective spin j, we invoke a general symmetry analysis, 
beginning with the non-Kramers case. Under time reversal, 
t z — > — T z jT x ' v — > t x ' v while j — > j. (Here, t z is by 
convention the component of r pseudospin along the local 
3-fold axis at the R-site.) This property leads to a remark- 
able simplification of the coupling - only t z couples to the Ir 
spin j. We consider the nearest-neighbor (NN) R-Ir exchange, 
which, due to space group symmetry, is parametrized by two 
couplings c\, c 2 . For the single Ir site labeled Irl in Fig.[T] the 
f-d exchange is 

n fd = [ Cl r| - c 2 (r| + rl)]jf + [ Cl r| - c 2 (r| + T*)]ft 
+ [cir| - c 2 (r z + ri)]jf + [2 O 2', 3 O 3', 4 o 4'], 

(1) 




FIG. 1. (Color online). A projective view of R2lr207 in the (111) 
plane. Left: the neighboring Ir (in dark red) and R (in light blue) 
tetrahedra. "1,2,3,4" label the four sublattices. Ir/R atoms are marked 
with big/small (red/blue) circles. Empty/Dark/Light circles indicate 
that the atoms are below/above/in the (111) plane. Right: an IrRg 
complex singled out from the left figure. 

In the Kramers case, the Ising part of the f-d exchange 
(coupling of t z to j) is identical. Transverse exchange involv- 
ing t x ' v is also permitted by time-reversal symmetry. Even 
in the non-Kramers case, while r x ' v does not couple to the Ir 
effective spin, it can couple to the Ir charge density. In both 
cases we ignore these transverse couplings, both for simplic- 
ity, and because they may be suppressed by strong easy-axis 
anisotropy of the /-moments along local [111] axes, which is 



present in many pyrochlore oxides. However, effects of trans- 
verse exchange may be important, and will be an interesting 
topic for future study. 

For the Ir subsystem, we follow Ref. 6 and include both the 
indirect hopping of 5d electrons through oxygen, and direct 
hopping between Ir sites, using the following Hubbard model, 

(rr') r 

where d\ a is the electron creation operator, with a I 
labeling the j z = 1/2,-1/2 states at site r, and n r ^ a — 
sum is over NN pairs of Ir sites. The direct 
hoppings (711/ ) involve two parameters|6), t„ and t n , that de- 
scribe the (7 and tt bonding, respectively. To be specific, we 
set = —\t a throughout the paper. The indirect hopping 
(Trfi) only has one hopping parameter which we denote as 

m. 

The R local moments can couple to each other either via 
superexchange through intermediate atoms, by dipole-dipole 
interaction, or by the RKKY exchange mediated by Ir elec- 
trons. Dipole-dipole interactions may play an important role 
for R (= Gd, Tb, Dy, Ho) where a large local magnetic moment 
is observed! 24]. RKKY exchange is likely to be the domi- 
nant exchange for the other compounds, as the Curie-Weiss 
temperatures in many of the isostructural insulating materials 
R2Sn2C»7 ||3D are of much lower magnitude than the corre- 
ponding iridates. For example, The Curie- Weiss temperatures 
are Q CW « -0.35K in Pr 2 Sn 2 7 and Q C w ~ -10KE1 
or -20KJ321 in Pr 2 Ir 2 7 . For the R = Nd compounds, 
®cw ~ -0.31K in the stannate J3l] and Q C w ~ -19K 
in the iridate ll24ll . 

From the above analysis, we obtain our minimal model 
for R 2 Ir 2 07, which includes the R-Ir exchange coupling in 
Eq. ([1} and the Ir-Ir hopping and interaction Eq. Q, 

rtmin = rlfd, + Hir- (3) 

To analyze the phase diagram, we start with the tight- 
binding model of the Ir subsystem. Following Ref. [6 1, a semi- 
metallic phase is obtained for — 1.674 <t a < — 0.67<. Other- 
wise, a strong topological band insulator (STI) with topologi- 
cal class (1;000) l33H35ll is obtained. In the semi-metal phase, 
at the r point there is a quadratic band touching (protected by 
cubic symmetry) at the Fermi energy (Ep). There are also 
non-dispersing bands at Ep along the T-L lines; this feature 
is a consequence of fine-tuning; it can be removed by adding a 
weak next-nearest-neighbor hopping (£') [6 |. The low-energy 
features of the band structure agree rather well with the first 
principles calculation for Y 2 Ir 2 07 [5 , 36 1, with the differences 
that the quadratic T-point touching is below Ep and the T-L 
lines have a small dispersion. 

We now add the f-d exchange to the Ir tight-binding model. 
Due to the Ising form of the coupling, the model does not con- 
tain quantum fluctuations of the /-moments, and reduces to a 
free fermion problem for any fixed configuration of localized 
moments. Finding the ground state amounts to finding the 
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minimum-energy configuration of local moments. Moreover, 
certainly c\, C2 <C t, so the f-d exchange can be treated per- 
turbatively, and the leading effect is to generate a RKKY ex- 
change between the /-moments. We find [30] that beyond 4th 
neighbors the RKKY exchange becomes significantly smaller, 
so we keep only up through 4th-neighbor exchange. Using the 
Luttinger-Tizsa method ll37ll . we find that the ground state of 
the truncated RKKY exchange has a q = magnetic order 
except in the light shaded regions of Fig. |2ja,fo). In the light 
shaded regions, the hard-spin constraint cannot be satisfied, 
and the nature of the ground state is not presently clear. How- 
ever, it is likely that the q = magnetic order extends at least 
somewhat into to the light shaded regions. 

Without losing any generality, we can simply focus on 
the case with c x > and define $ = tan _1 (c2/ci) and 
c = y/ci + cf. As shown in Fig.^a, b), for most of parame- 
ter space, "all-in all-out" magnetic order is favored, where ev- 
ery tetrahedron of neighboring R sites has either all t z point- 
ing in {i.e. toward the tetrahedron center), or all pointing out. 
In the dark shaded region, q = "two-in two-out" magnetic 
order is obtained, where on every tetrahedron, two t z point in 
and two point out. [The q = two-in two-out state also has 
lowest energy, at least among q = states, in the vertically 
hatched regions.] Since no ferromagnetic state is observed in 
any R2lr2C>7, we restrict our discussion to all-in all-out mag- 
netic order. Such order of the R subsystem also induces all- 
in all-out magnetic order in the Ir subsystem via the f-d ex- 
change, which acts as a local magnetic field modifying the Ir 
band structure. 

As shown in Fig. |2fa,&), with an f-d exchange, the STI 
phase is converted into an AI for small c/t and then into a 
WSM. The semi-metal phase immediately becomes a WSM. 
In this case, it is necessary to add a very small t' to stabilize 
the WSM 0; without it the flat T-L lines of the semi -metal 
remain at Ep. 

Returning to the STI case, when time-reversal symmetry is 
broken but inversion symmetry is preserved, as for the all-in 
all-out state, the magnetoelectric response parameter 9 is still 
quantized. We can express 9 in terms of the number n (k) of 
filled odd parity states at the time reversal invariant momenta 
(TRIMs) MM, 

- = l Y, »o(k) (mod 2). (4) 

ke TRIMs 

As shown in Fig. |2|a), a region of AI phase with 9 = ir is 
obtained in the vicinity of the STI phase. The band structures 
of the STI and the corresponding AI induced by f-d exchange 
are shown in Fig.|5Ja, b). 

The Fermi surface of WSM is composed of the Weyl points, 
around which the spectrum is linear and gapless. The effec- 
tive Hamiltonian near a Weyl point is obtained by expanding 
the full Hamiltonian at the Weyl point with only two relevant 
bands involved[5|, 

"Hwsm = E F + v • q + ^(v, • q)cr i7 (5) 




FIG. 2. (Color online). Phase diagrams of the STI with t a = — 2t 
(in a) and the semi-metallic phase with t a = —t (in b) after in- 
cluding f-d exchange. See text for further discussion of the phase 
diagrams. In (c, d), the Curie-Weiss temperature Ocw (solid line) 
and the mean-field ordering temperature T c (dashed/dotted line) of 
the truncated exchange model (up to 4th neighbor for R system) are 
plotted against $. T c is obtained by restricting to the q = mag- 
netic ordering|30|. T c for the all-in all-out (two-in two-out) state is 
the dashed (dotted) curve. In (e), we plot the c%-t a phase diagram for 
Ir subsystem where the magnetic order is all-in all-out. In the figure, 
NI = normal insulator. 

where q is measured relative to the Weyl point, and <jj are 
three Pauli matrices. A chirality C, which is defined as 
C = sgn\y\ ■ (V2 X V3)], can be assigned to each Weyl point. 
Weyl points can only be created or annihilated in pairs of the 
opposite chirality. We locate the Weyl points explicitly by ex- 
amining the spectrum and density of states, which shows the 
characteristic (E — Ep) 2 scaling. As c increases, the 8 Weyl 
points of the WSM in Fig.|2j6) are created at T and annihilated 
at the L points, and lie along the T-L lines [Fig.[3jd)]. 

The Weyl points are sources and sinks of the Berry curva- 
ture. As argued in Ref. Weyl points with opposite chiral- 
ities should be connected by gapless surface states that form 
open Fermi arcs in the surface Brillouin zone (SBZ). As the 
Weyl points in this WSM phase appear on a line with high 
symmetries, the Weyl points of opposite chiralities project to 
the same point on the SBZ for (001) surfaces. Therefore, no 
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FIG. 3. (Color online). The evolution of the Ir electron band struc- 
ture. The energy unit is set to t. The dashed (red) line is Ef- (a) is 
STI phase with t a = — 2t and ci = C2 = 0. (b) is AI phase with 
t a = —It and ci — 0,C2 = 0.05i. (c) is semi-metallic phase with 
t a = —t and ci = C2 = 0. (d) is WSM phase with t a = — t and 
ci = 0,c 2 = 0.05tf38). 



surface states are expected for these surfaces. For (111) and 
(110) surfaces, however, Weyl points with opposite chiralities 
project to different points in SBZ, leading to gapless Fermi 
arcs on the surface. 

We now include the correlation on the Ir subsystem. With- 
out the f-d exchange, the Hartree mean-field analysis of the 
extended Hubbard model on Ir subsystem gives a rather nar- 
row region of WSM phase [6'|. With f-d exchange, the RKKY 
exchange induces magnetic order even with weak correlations 
and thus gives a WSM phase. We start from the metallic phase 
and include both the Hubbard interaction (via the same mean- 
field approach used in Ref. [6 1) and f-d exchange. Since both 
correlation and RKKY interaction alone give an all-in all-out 
magnetic order on Ir system, we restrict attention to the all-in 
all-out state. We decouple the Hubbard interaction, 
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FIG. 4. (Color online). The mean field phase diagram at t a = —t 
and ci = 0. The (red) thick line is the narrow WSM phase induced 
by Ir correlation. The dashed line is the nonmagnetic Ir semi-metal 
phase. The (green) region is the broad WSM phase after including 
f-d exchange. All-in all-out magnetic order is present throughout 
the phase diagram. 



A much-studied compound in the R2lr2C>7 series is R = Pr, 
which exhibits a metallic ground state and anomalous Hall ef- 
fect without any observable magnetic order for 0.3 K < T < 
1.5 K(20l[32]|3l). The Pr local moments freeze at T < 0.3 K. 
Our theoretical analysis focuses on the all-in all-out region of 
the phase diagrams in Fig. [2] Pr2li"207 may be located in or 
near the light shaded region of Fig. (2^6), where the RKKY 
exchange may be more frustrated. Moreover, quantum fluctu- 
ations may be important for R = Pr, perhaps originating from 
the transverse part of the f-d exchange, from superexchange 
between the /-moments [14|, or both. Such quantum fluctu- 
ations may suppress the magnetic order, perhaps cooperating 
with classical frustration. 



Nd2lr207 has been observed to order magnetically below 
the metal-insulator transition at T c = 36 K EW3\ [291 . Re- 
cent neutron scattering experiments [23] on Nd 2 Ir 2 07 suggest 
an all-in all-out spin configuration for both Nd and Ir systems. 
The data was interpreted in terms of a static local field set- 
ting in for the Nd moments immediately upon cooling through 
T c , while Nd magnetic order is not detected until 15 K. An- 
other experiment! 22 1 finds that applying pressure to Nd 2 Ir 2 07 
suppresses the metal-insulator transition and leads to a metal- 
lic state that exhibits a negative magnetoresistance. Since 
9cw = — 19 K in Nd2li"207, according to Fig. (2jc,<i) we 
expect Nd2lr2©7 to be in the region with $ < in order to 
have an antiferromagnetic &cw- Moreover, T c > \®cw\ an d 
the all-in all-out ordering are also consistent with Fig.|2];c,(i). 
The metal-insulator transition and magnetic order are likely 
driven by a combination of RKKY exchange and Ir magnetic 
exchange. Indeed, Ref. f23l interpreted their results in terms 
of Ir magnetic order setting in at T c , with Nd order setting 
in only at 15 K (driven by f-d exchange) due to the presum- 
ably small Ir ordered moment. It would be interesting to see 
whether our model can indeed produce such behavior at fi- 
nite temperature. Assuming all-in all-out ordering is indeed 
present, the Ir electrons may realize a WSM phase within the 
magnetically ordered state. Indeed, the transport properties of 
the WSM are close to those of an insulator[40], and appear 
consistent with the resistivity of the Nd compound below T c . 
Applying pressure to Nd2lr2C>7 [22] may modify the Ir elec- 
tron hopping, generating a metallic state. 

To summarize, we have constructed and analyzed a mini- 
mal model accounting for f-d exchange in rare-earth-based 
pyrochlore iridates R 2 Ir 2 07. With this model we find broad 
regions of Axion insulator and Weyl semi-metal phases. 
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SUPPLEMENTARY MATERIAL 



f-d exchange between the R moment and Ir moment 

Both Ir and R sublattices are pyrochlore lattices, and can be viewed as FCC lattices with a four-site basis. We use the BCC 
primitive vectors 

b x = (0,1/2, 1/2), (7) 
b 2 = (1/2,0,1/2), (8) 
b 3 = (1/2, 1/2,0). (9) 

For the Ir pyrochlore lattice, we choose the following basis vectors, 



in = 1/4(0,0,0), (10) 

Ir 2 = l/4(0,1,1), (11) 

Ir 3 = 1/4(1,0,1), (12) 

Ir 4 = 1/4(1, 1,0). (13) 
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For the R pyrochlore lattice, we choose 

Ri = (0,1/2,0), (14) 

R 2 = (0,3/4,1/4), (15) 

R 3 = (1/4, 1/2, 1/4), (16) 

R 4 = (1/4,3/4,0). (17) 
The symmetry-allowed f-d exchange between the R moment and the Ir moment can be written as 

Ufd = E T ^ r ) h [32 ( r ) + J2 (r b 2 + b 3 ) + if (r b 2 + b 3 ) + if (r + bi - b 2 ) + if (r) 

r 

+ if (r + bi - b 2 )] + c 2 [jf (r) + if (r - b 2 + b 3 ) + jl(r) + if (r - b 2 + b 3 ) + if (r - b 2 + b 3 ) 
+ if (r + bi - b 2 ) + if (r b 2 + b 3 ) + if (r + bi - b 2 ) + if (r) + if (r + bi - b 2 ) + if (r) 

+ jj(r + bi-b 2 )]" 

+ 7- 2 *(r) [ci [if (r + b x ) + if (r + b : b 2 + b 3 ) - if (r + bi) - if (r + b x b 2 ) 

- if (r + bi - b 2 ) - if (r + ^ b 2 + b 3 )] + ca [ - if (r + bi) - if (r + bi - b 2 + b 3 ) - if (r + b^ 

- if (r + b! b 2 + b 3 ) + if (r + bi - b 2 ) + if (r + bi) - if (r + bi - b 2 ) - if (r + bi) 

+ if (r + bi - b 2 ) + jf (r + bi - b 2 + b 3 ) - if (r + bi - b 2 ) - if (r + b 1 -b 2 + b 3 )] 

+ rf (r) [ Cl [if (r + b 3 ) + if (r + bi) - if (r) - if (r + b 3 ) - if ;G (r) - if (r + b,)] 
+ ca [ - if (r + b 3 ) - if (r + bx) - if (r + b 3 ) - if (r + bi) - if (r) - if (r + b 3 ) 
+ if (r) + i|(b 3 ) + if (r) + if (r + b x ) - if (r) - if (r + br)]" 
+ rf [ci [if (r + b 3 ) + if (r + bx - b 2 + b 3 ) - if (r - b 2 + b 3 ) - if (r - b 2 + b 3 ) 

- if (r b 2 + b 3 ) - if (r + b! - b 2 + b 3 )] + en [ - if (r + b 3 ) - if (r + bi b 2 + b 3 ) - if (r + b 3 ) 

- if (r + bx - b 2 + b 3 ) - if (r - b 2 + b 3 ) - if (r + b 3 ) + if (r - b 2 + b 3 ) + if (r + b 3 ) 

- if (r b 2 + b 3 ) - if (r + bx - b 2 + b 3 ) + if (r - b 2 + b 3 ) + if (r + b x - b 2 + b 3 )] ] . (18) 

Here rf (r) is the pseudospin for the R site in unit cell labeled by r with sublattice index i = 1, . . . , 4. (Recall that t z is the 
component of pseudospin along the local 3-fold axis at the R site.) Similarly, if (r) is the ji = x,y 7 z component of the Ir 
effective spin (in global cubic coordinates), for the Ir site in unit cell r at sublattice i. c\ and c 2 are the two parameters allowed 
by space group symmetries. 

It is convenient below to work in a local coordinate system for the Ir effective moment. So, we transform the Ir effective spin 
from glocal to local Ir0 6 octahedral coordinate system by 



if(r)=^ri- L (r), 



(19) 



where is the Ir effective spin in the local IrC>6 octahedral coordinate system. The transformation matrices for four 

sublattices are 



K x = 



n 2 = 



n 3 



TZ 4 



2/3 


-1/3 -2/3 


-1/3 


2/3 -2/3 


2/3 


2/3 1/3 


2/3 


2/3 1/3" 


-2/3 


1/3 2/3 


1/3 


-2/3 2/3 


1/3 


-2/3 2/3" 


2/3 


2/3 1/3 


-2/3 


1/3 2/3 


1/3 


-2/3 2/3 


-2/3 


-2/3 -1/3 


2/3 


-1/3 -2/3 



(20) 



(21) 



(22) 



(23) 
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RKKY interaction between the R moments 

We consider the Hamiltonian 

H=n b + H fd , (24) 
where T~Li r is the tight-binding Hamiltonian without the Hubbard interaction, 

= ^ (7^r',o!/3 + 7?r',a0)diadr'l3- (25) 
(rr') 

The f-d exchange H / d can be written compactly as 

n f« = E E - r ')ijV)- (26) 

r,r' 

Here, and below, sums on repeated indices are implied. The exchange is characterized by the function (r), which can be read 
off from Eq. ( fT8| ). 

Diagonalizing the tight-binding Hamiltonian T~Li r gives 



Hlr = ^ £ kAc| (:A C k A, (27) 



k,A 

where A labels the eight eigenstates for each k, and 

4 M -E c kA^A, M (k). (28) 

A 

We now treat Hfd as a perturbation to the Ir tight-binding Hamiltonian W^. By standard second-order perturbation theory, we 
obtain the effective RKKY exchange between the local moments of the R subsystem, 

n RKKY - ^ Eo^e, ' (29) 

i 

in which, |$ ) an d Eo is the ground state and the corresponding energy, and and E{ is the ith excited state and the 
corresponding energy. This result can be written 

n RK KY = \ J2 J M 2 ( r i- r 2K(riK(r 2 ) (30) 

riii,r 2 l2 

with J ili2 (n - r 2 ) = J ni2 (ri - r 2 ) + Ji 2 i x (r 2 - ri), where 

Jixi 2 (ri - r 2 ) - 1 2^ ■'nii^i ~ r iJ/i 2 j 2 t r 2 - raJ/c^ tt ja cr Qi/3i cr a2/32 2^ e 

r^,r 2 c ki,k 2 

MAi jiai (Ki )W A • ^ (ki )^a 2 , j 2 a 2 (k 2 )W A2 Jlft (k 2 J — — . (31) 

e \i\ k l) ~ £ A 2 (K 2 j 

Here, u v is the Pauli matrix, N c is the number of unit cell, and Q(x) the Heaviside step function. 

Since the RKKY exchanges are quadratic in c\ and Co, we can simply focus on the regime with c\ > and define an angle 
$ = tan _1 (c 2 /ci) and c = y/ cf + cf. As shown in Fig.pl we compute the RKKY exchanges of a finite system with 10 x 10 x 10 
unit cells for both the STI phase with t a = —2t and the semi-metal phase with t a = —t. Clearly, the RKKY exchanges beyond 
4th neighbor become significantly smaller. We now focus on the truncated RKKY exchanges upto 4th neighbor. The Curie-Weiss 
temperature is found to be, 



Ocw = M+2J2-3J 3 -34 + 2Jt). (32) 



The transition from the high temperature paramagnetic phase to low temperature ordered phase is continuous at mean field 
level. The transition from paramagnetic phase to all-in all-out state is of three dimensional Ising type. We can then determine 
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the critical temperatures for this transition by the usual condition of marginal stability (vanishing of the quadratic term in the 
Landau theory) of the free energy. Restricting to q = magnetic order, we find that the critical temperatures for all-in all-out 
and two-in two-out states are given by 

T c (all-in all-out) = -(-J 1 - 2J 2 -J3-J3- 2J 4 ), (33) 

T c (two-in two-out) = -(Jx + 2J 2 - 3J 3 - 3Jg + 2 J 4 ), (34) 

respectively. The actual transition temperature is determined by the higher one of the above transition temperatures. 



^/\c 2 /t\ Ji/\c 2 /t\ 




FIG. 5. (Color online) The dependence of RKKY exchanges on the angle $. The results are obtained for a finite system with 10 x 10 x 10 
unit cells. Ji is the ith neighbor exchange. J3 is the exchange on 3rd neighbor bonds that have one lattice site at the midpoint of the bonds. 
J3 is the exchange on 3rd neighbor bonds that have no lattice site at the midpoint of the bonds. J7 is the exchange on 7th neighbor bonds that 
have two additional lattice sites along the bonds. Jj is the exchange on 7th neighbor bonds that have no lattice sites along the bonds. Upper: 
STI phase with t a — —2t. Lower: semi-metal phase with t a — —t. 



